


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1973 


Design on the frequency parameter plane 
with application to roll stabilization of ship. 


Vachanaratana, Manop. 


Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/16559 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
) (Ч D U DLE Y research materials and institutional publications created by the NPS community. 
4 : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


ІШ KNOX appointed — and published — scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 











AVAL POSTGRADUATE SCHOOD 
MON САТТЕ, 92940 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


Design on the Frequency Parameter Plane 
with 
Application to Roll Stabilization of Ship 


by 


Manop Vachanaratana 


Thesis Advisor: G. J. Thaler | 


December 1973 


Approved for public release; distribution unlincted. 








Design on the Frequency Parameter Plane 
with 
Application to Roll Stabilization of Ship 


by 


Manop Vachanaratana 
Lieutenant Juniór Grade, Royal Thai Navy 
B.S., Royal Thai Naval Academy, 1964 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1973 





Library 
Naval Postrraduate School 
Monterey, California 939^ 


ABSTRACT 


Frequency parameter plane techniques are developed for 
the analysis and synthesis of both linear and nonlinear 
systems. The transmission function of a dynamic system 
contains two parameters which are adjusted to provide 
desired performance. These parameters represent elements 
of the system. 

The four dimensional paremeter space relating magnitude, 
frequency, and the two parameters may be represented by seven 
two dimensional projections of the parametric curves. 
Computer programs are developed for the computation, and 
graphs of results are obtained. Interpretations of the graphs 


can be used for design of systems such as tank stabilizers. 
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I. INTRODUCTION 


The analysis and synthesis of control systems may be 
accomplished using the frequency parameter plane technique. 
A primary consideration in control system design is the 
adjustment of the system frequency response. To achieve 
this, the frequency parameter plane technique may be employed 
and is based on a graphical procedure for considering four 
parameters, alpha, beta, frequency and magnitude of the 
transmission functions. The objective is to choose values 
for alpha and beta (or functional relationships if the design 
is to be nonlinear) to provide the desired frequency response 
characteristics. This choice is possible through proper use 
of the frequency parameter plane curve families. 

The goal of this thesis is design using analysis techniques. 
The parameter frequency response method is used to synthesis 
and design systems. A steam-engine governor is used in the 
analysis as an illustration. The antiroll tank stabilizer of 
a ship is also designed using the frequency parameter plane, 
for both passive and active tanks. The tank size and active 
device characteristics are represented by the parameters, 
which are varied. Parameter plane curve families are used 


in the design procedure presented in later chapters. 
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II. PARAMETER PLANE ANALYSIS 


Dynamic systems may be either linear or nonlinear. 
Parameter plane analysis permits us to study dynamic behavior 
by treating system parameters as variables. One first 
determines the characteristic polynomial of the system which 


may be in the form 


k N 


5 = ME ta S rac ой: t a 


N 
Е(5) з 7 а, 1-1 


К-0 
(2.1) 


where ay = coefficients which are functions of the parameters. 


Е о а - 1 = 39 
5 tjo 20,330, ut ве 


By substituting the value of S into equation 2.1, 


equation 2.2 is obtained. 


N : 
Е (5) + 7 а, m e JK8 (2.2) 


Since equation 2.2 has complex coefficients it can be converted 
into two equations by requiring reals and imaginaries to be 


zero independently. 


Note that S* - р jiko | ық (Cos ke + j Sin ke) 


1 
Е 
(D 

| 


ша 
where X, - wrCos КӨ = вСо (К Со5710) - (-1) окт: (с) 


k 26811 kg = oKS i n(k С05712) = (-1)Kwku, (2) 


< 
li 
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Тұ(с) and 0,(5) аге Сперузлеу polynomials with values 
as shown in Table 1, 2, and 3. 


It follows that 


a, EC n (1 = 0 


1 Mm 


k=0 


ap[(-1) “ou, (g)] = 0 (2.3) 


ж” 
п мМ 2 
© 


are simultaneous equations which may be solved for two 
undefined parameters in the а“ 
For many practica] problems the algebraic form of the 


coefficients is one of the following two cases: 


b, a + с В т 4, 
= bra + c, B т n t 4, 


Case 1 can be solved by Cramer's rule and the results 


for a and 8 are 


Е 09 
T BoD, ~ D5B4 
21% - BoC, 
where 
N N 
а К К . E k k . 
By = Е (1) Tele): Ba Ep (1) CB 0,0) 
N N 
К 
С: -(20(-1) сү, T (5) ; С» 2 Ч (-1) Кско Ки (с) 


N k К Ч 
р, m dr Т (е); ре ce 


k=0 
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Case 2 by substituting equation 2.3 rearranges into 


the form. 
aB, + BC, + aBH, + D, = 0 ШКО ОООО у... (2.4) 


where the definitions of case ] still apply, and 
N N 
k К К К 
= > (1) ют (5): Н. мо (0) h,o, U,(Cz) 
1 k=0 kn k 2 k=0 kn “k 


Eliminating 8 the solutions are obtained 


. -etye?- Fac 
91 „2 ш 2а € & ә ә ® à ® à * * 9 * * * * е (2.5a) 
and corresponding values for ß are 
ME. AZ um 
1:52 Но 2С1 Но 276, 


Eliminating a, the following solution for 8 


вур = 06 ................ (2.6а) 
апа the corresponding values for o are 
2-6 30 USB *D 
а = Д = hii 222211. (2.6b) 
к 1.21 25172-222 


where the notation is the following 


а = В.Н, S с = CD» - CoD, 
b= CoH, - СІН; d = B,D, - BoD, 
e = ВС, - B1C5*H4 D-H5D, 
Ши 12:12 17221727) 
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As can be seen, there are generally two pairs of solu- 
tions (a,8) in terms of ш, апа t. The solutions are given 
by equations 2.5a and 2.5b provided a # 0. If a = Тоог 
solutions are given by equations 2.6a and 2.6b provided 
b #0. When a = b = 0 case 2 reduces to case 1. 


Illustrations I. A feedback control system is as shown 


=), 






The closed-loop transfer function is: 


K 
T(s) - | 


53 + 352 + (24K K5)S * K 


2) 1 


The characteristic equation of the system is 


E 35° + (2*K.K 5-4 K4 = Q 


2) 


! 
> 
> 


Let а 
В = К 


The equation becomes: 


3 


шэг 


From the parameter plane method 
В. =- u. В, = -wn 


С. = 1 СЕО 
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(483-35) (-ш 3) + (252-1) (3ш 2) + (-2tu,) 


O 
i 


D, - (422-1) (-4,?) + (25) (34,2) + (-2u,) 


50, the solutions for a and В аге 


а --4u * c + бш 5 t ше - 


sz Зиг? Е 2сы З 


The z-curves and and w -curves on the a-8 plane are shown 
on Figure 2.1. 
Illustration II. The system is the same as illustration I 


let К, = В 


Ko а 
The characteristic equation becomes: 


со + 35° + (2+а8)5 + B = 0 


Then 

В, - р B, = 0 

Ci - 1 C5 = 0 

Hy = -to, Н, iz 
u 3 3 2 2 

р, = (35-45) un + (67 -3) u, ^ 2 Zu, 
s 2 3 2 

D, = (1-47) шд + бсш - 20, 


And the other definitions are 


a - 0; b= ws с = 0; 
2-2 2 


d- -Au Pt + 6 w 2 


2 3 
n^ t (30, -20,) 2 - Зи" 


ХЭГ 721777 3,7 
n n 
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So only one solution for a and 8 in equations 2.6a and 2.6b 


is obtained: 


= 2 З 
В = 30, - гш С 
B 222 2 2 
a = 40,5 - 6w,5 - on өл 2 _ 25 Е Qu 7 2 
3 2 0) В 
20.5 - 3o, n 


The curves are shown on Figure 2.2 
Consider the curves on Figures 2.1 and 2.2: the roots of 
the characteristic equation are defined as functions of a and 


8. If we choose g = 40 and a = 20 on Figure 2.1, then 


4.5 


0) 
n 


С o Sp 


therefore, the roots of the system are: 
5 = - 0.5 + j4.47 

and the corresponding value for the real root is -2. In like 
manner on Figure 2.2, ¡if we desire roots at С = 0.111 and 
харе 7.5 we found a = 0.5 апа р - 40. 

On the other hand, if z = 0.2 and б. 5 the roots are 
-1:j4.9 and -l. From the Figure we find: 

B = К, = 25 
25 (Figure 2.1) or 1.0 (Figure 2.2) therefore K 


а = 1.0 


2 
Parameter plane methods are very convenient. The а-8 
plane can be used to find the system parameters, and/or roots 
of the characteristic equation at any defined values of the 
parameter. The stability of the system can be investigated 


at the line c=0 which corresponds to the imaginary axis on 


the S-plane. 
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Figure 2.1 Root Parameter Plane Curves for illustration 1. 
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III. PARAMETER FREQUENCY RESPONSE TECHNIQUES 


A. INTRODUCTION 

Both linear and non linear dynamic systems may be sub- 
jected to periodic inputs. It is important therefore that 
we be able to analyze such systems in the frequency domain. 
If the system is linear there also exists a correlation 
between frequency response and time response as guaranteed 
by the Fourier integral which permits us to estimate the 
transient response from the frequency response. Thus the 
frequency domain may be used to design systems for non- 
periodic inputs. In this chapter frequency response analysis 
of linear systems on a frequency parameter plane is discussed 


in detail. 


B. ANALYSIS 
A linear system can be described mathematically by a 


transmission function: 


Ms) = 909) 


where the system may be either open-loop or closed-loop and 
N(s) and D(s) are polynomials. In polynomial notation 


а 58 
=0 
Т(5) = 


Ь $9 
-0 4 


о мМ зія м = 


where, a's and b's are coefficients which are functions of 


the system parameters. 
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If we choose two systems parameters a and B and if we 
consider only sinusoidal signals the frequency response 


transmission function becomes 


5 2 N(jw,a,ß) 
TK jw) = D( jo, 0,8) 
я ag * ajo t а „(јш) 2+ ДЕЛ р + а (Зо)" 
by + byjw * bo(Ju) ^... sese * b (ju) " 
: 3 
_Lag - а уш + адо ар ар ти 153|Гаүо-а, Pag... 
[bo - b u^ + byw ОЪ ]+3[b,0-b2u°+b.u°-... 
In abbreviated notation 
М. (ш,а,В) + jN (w,a,8) 
T(jw) = ое 90205000 
e 3 э J 0 9 9 
Similarly 
И с, № (ш,а,В) - 3 N (w,a,ß) 
-497 7 D. (u,a,B) - 10 (4.88) 
Їр о „а „В 2 Bo Gc sp 
Let the magnitude be represented by M: 
М = |Т(3ш)| 
Cy ; a : : 
шиг OL 70)! 
therefore 
> IN, (w.a,8)]* - [N (wsa,8) 1° 
vts 2 2 
ГО (юьо, В) 1 т [D,¿(0,0,8)] 
from which 
2 4 2 3 2 
E f [ар-ар ғадш -...... 1 + Гал ш-азо Боолын 1 
[b,-b we +b vm ie + [b,w-b 3+ 1° 
0 2 4 аа Га ә s 3 вас оао 
2m 820,0 2п-1 к, ү 
L (21) 5 ( -1) а, в 
_ k=0 (even) k=] (odd) 
| On q д1- | 4-1 
У AA MOL 1 
q=0 (even) “q q=1 (odd) q 





Then 


where, the ms and bq 5 are functions of a and В. So, there 
are four variables M, o, a, and 8 in this equation. The 
equation can be plotted by fixing one variable and using one 
variable as a parameter. Seven different combinations are 
possible. 


Example. Consider the system of a Steam engine Governor 


[9] for which the transmission function can be written as: 





СС, Р 
$ - Im 
$c 54 (Са+ Се) 53 + (Ке+ Ка+ СеСа) 5 2, (те e qe ид % 2)58-6 M 
m I Itam lm lan и тип m 


where @ = instantaneous speed response. 
$c - command speed. 


m = equivalent mass of the governor sleeve. 


Eg - damping coefficient at the governor sleeve. 
K = stiffness of governor spring. 
С, - increase in the upward force on the governor 


sleeve. 
с, = increase in steam torque. 
I = equivalent moment of inertia of the engine. 


С. - damping coefficient of the engine. 
K 


И 


stiffness of engine spring. 
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From studying the mechanical characteristics [10], 
possible values of the parameters of the system are sub- 
stituted in the equation. Let Ka and с, be varied by adjusting 
the governor spring and the valves for steam flow respectively. 


Then the equation becomes: 


я 626«5 
[» 3 
фо 54.10.553+(30+258)52+(137.58+6260+12.5)5+62.58 
where а = С) 
ми 
р 


then 


м2 » ! [62600] 
by 2 
) 


| «"-(308258)ш5562.581541-10.50"-(137.5846260412.5)01 

From the curves of Figure 3.1 several useful bits of 
information can be obtained by inspection: 

(a) At any intersection between two constant w curves 
the coordinates define the a-8 pair, and the values of w of 
the intersecting curves define the bandwidth (since the mag- 
nitude is -3 db). 

(b) A locus of constant w define a set of a, B pairs 
which are potentially useful in self adaptive applications, 
i.e., if, in a given system, 8 is subject to change, values 
of a can be determined such that a magnitude M = М. can be 
maintained at w = w}, by determining the values assumed by 8 
and providing a suitable adjustment in a. 

(c) If a number of families of curves are generated with 
coordinates as in Figure 3.1 but with different M values, the 
entire frequency response can be determined for any chosen 


ieee pair. 
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From the curves of Figure 3.2a, b and c, it can be seen 
that: 


(a) Each frequency has only one a, ß pair which gives 


the maximum magnitude. 


i.e., atw = 10, а = 1.0, 8 = 3.0 (Figure 3.2a) 


12, а 1.4, 8 4.5 (Figure 3.2b) 


u) 


о = 14, a= 1.75, В 6.6 (Figure 3.2c) 

(b) The locus of the maximum magnitude point is a 
Straight line when frequency is increased. 

(c) The axis which is drawn through the ellipses gives 
the values of 8 (and/or a) and magnitude which are coordinates 
of points on the envelope of the a versus magnitude curve as 
in Figure 3.3. 

On Figure 3.3 are curves of o versus magnitude. These 
curves as drawn show the envelope of the family of curves 
that corresponds to the axis of the ellipses on Figure 3.2. 
This envelope bounds all possible magnitudes at a frequency 
of 12 rad/sec, for all 8. 

Figure 3.4, which presents a family of constant magnitude 
curves on a versus frequency coordinates (8 = constant) pro- 
vides a global view of the possible frequency responses of 
the system: 

(a) We can plot the Bode diagram by reading the frequency 
and magnitude at the intersections along any line of constant a. 

(b) The maximum magnitude for the given values of 8 is 
found at the center of the family of ellipses, and also 


increases as a straight line when 8 is increased. 
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(c) At constant 8 (855.0), the curve of magni tude 
is -3 db., it shows the bandwidth is between the intersection 
of the curve and any line of constant a, such as ata = 2 
the bandwidth is 10.5 rad/sec. 

A plot of 8 versus frequency with constant a may be used 
in the same way. 

Curves of a versus magnitude are shown on Figure 3.5 
with 8 constant and frequency the family parameter. 

(a) The Bode diagram can be plotted at any a by evaluating 
magnitude and frequency along a horizontal line. 

(b) At the intersections of two frequency curves the 
Bode diagram for that value of a will have the same magnitude 
at each of the frequencies. 
(c) Each frequency curve defines one value of a (B is 
constant) for which the magnitude is maximum at that frequency. 
(d) If two frequency curves intersect each other at 
the magnitude -3 db., the bandwidth at that value of a is 
defined by the two frequencies. 

Figure 3.6 plots a versus frequency with g as the para- 
meter. Each ß-curve is intersected twice by each a line, 
thus defining two values of frequency for the one value of 
the a, 8 pair, these two values of frequency give the band- 
width of the system (since the magnitude is -3 db). 

Figure 3.7 shows a family of Bode diagrams with a as the 
family parameter and 8 constant. It may be used to select a 
value of a, or to define a desired variation in a as a function 


of frequency. 
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Each curve family is seen to contain useful information. 
Depnding on the type of information needed some curve families 


may be more useful than others. 


C: STABILITY ON THE FREQUENCY PARAMETER PLANE 
Figure 3.8 shows the stability on the root parameter plane 


of the above system. The curve of єк = 0 is obtained from: 


, 2 1.248 - 0.49020 * 5.5410? « 0.199684" 
(ОБЩО = 62.5) 
g. We 304 
250^ - 62.5 


In the typical system, the region of stability on the 
frequency parameter plane can be determined by using the 
above relations. 

The stability boundary (£20) curve has been added on 
ura 3.1, through Figure 3.6 and the shaded areas are the 
stable region. Synthesis and design with the frequency para- 
meter plane are considered only for stable systems. 

Each Figure has a different interpretation. Figure 3.1 
plots a versus 8 with magnitude constant at -3 db, and frequency 
as the family parameter. By choosing an operating point in the 
stable region at the intersection of two frequency curves, one 
guarantees both stability and bandwidth. Figure 3.2 uses 
frequency as the constant parameter and magnitude as the 
family parameter. This permits optimum choice of a and 8 


for desire magnitude and frequency. 
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The other curves can be considered in the same way. At 
any point on the Figures in the stable region, the parameters 
at that point can use to design the system and obtain values 


of magnitude and frequency as desired. 
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Figure 3.2a Frequency Parameter Plane with Stability Limits. 
Plot of a vs 8 with o constant at 10 rad/sec 
and M as family parameter. 
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Figure 3.2b Frequency Parameter Plane with Stability Limits. 


Plot of a vs 8 with w constant at 12 rad/sec 
and M as family parameter. 
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Figure 3.2c Frequency Parameter Plane with Stability Limits. 
Plot of a vs 8 with w constant at 14 rad/sec 
and M as family parameter. 
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Figure 3.4a Frequency Parameter Plane with Stability Limits. 
Plot of a vs w with B constant at 8-4 and M as 
family parameter. 
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Figure 3.4b Frequency Parameter Plane with Stability Limits. 
Plot of a vs w with 8 constant at 8-5 and M 
as family parameter 
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Figure 3.5 Frequency Parameter Plane with Stability Limits. 
Plot of a vs M with 8 constant at 8=5 at w as 
family parameter. 
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Figure 3.8 Root Parameter Plane Plot. 












IV. NONLINEAR SYSTEM TECHNIQUES 


Generally, nonlinear systems consist of linear parts 
which are described by linear differential equations and 
nonlinear parts which are described by nonlinear functions. 
If the system is nonlinear with an applied sinusoidal input, 
the signal at various points in the system may or may not be 
sinusoidal. In general the signals will be periodic, but 
may contain sub-harmonics and/or harmonics. 

The most widely used method in the stability analysis 
of nonlinear control systems linearized by the describing 
function technique is the Nyquist diagram. With the concept 
of the variable critical point, which depends upon the frequency 
and signal level at the input of nonlinearity, the Nyquist 
diagram permits a simple stability analysis ‘and investigation 
of sustained oscillation. 

However, the Nyquist diagram is inconvenient to apply to 
multiloop control systems. Difficulties also arise with the 
adjustment of the system parameter, since tlhe diagram should 
be plotted each time a parameter iS changed. Thus, this 
method can hardly be applied to design problems in which the 
describing function does not appear as a variable gain factor. 

Other methods may be used, namely the Hurwitz criterion, 
Root-locus techniques, Bode diagrams, Nichols charts, 
Mitrovic's method, approximation method, Popov method, and 
phase plane techniques, etc. which can be used to investigate 


the nonlinearity of the system. However the parameter plane 
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method is more convenient since the nonlinearity values can 
be represented by the parameters. 
Consider the sinusoidal system; 


Т(3ы) = №6058 


Then 
= 1 = N(jw,a,ß) 
а ITO) | ова 
Define one parameter as 8 to be constant and the other para- 


meter a represents the value of the nonlinear element. 50, 


н = [06825 


The equation has three variables: magnitude, frequency, and 





the nonlinear parameter o, which the parameter plane method 
is to plot. 
The block diagram of the nonlinear system may be the 


following: 


о G(s) 


% is the nonlinear element. 


G(s) represented the linear system. 
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Nonlinear elements such as square law devices, coulomb 
friction, dead space, saturation, backlash, relays, etc. may 
be represented by their describing function. The parameter 
is then used to represent any adjustable parameter in the 
describing function, and the resulting plots on the frequency 
parameter plane indicate the effect of the nonlinearity on 


the system behavior. 
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V. DYNAMIC SYSTEM DESIGN 


A. INTRODUCTION 

In the design of compensators for linear or nonlinear 
systems, linear system design to make the system stable uses 
the root locus, Bede plot, Nyquist and Nichol chart. However 
the parameter plane method is more convenient for both linear 
and nonlinear systems. In classical methods of the root 
locus and frequency domain (Bode diagram), the compensators 
are introduced to reshape the corresponding graphs of the 
original system in order to achieve desired system response 
characteristics. 

In the parameter plane approach the dtrect effects of 
the additional blocks and their connection is to change the 
characteristic equation, so that the additional adjustable 
parameters are conveniently distributed among the coefficients 
of the characteristic equation. This enables us to have 
control over the characteristic roots by varying the introduced 
Parameters. 

The parameter plane is used here to investigate and 
design the system for antiroll stabilization of a ship by 
means of passive tanks. The size and properties of the tank 
are designed. The problem is extended by including the pump 
System thus converting the system to an active stabilization 


system. 
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В. THE DYNAMICS OF A SHIP IN ROLL MODE 
A mathematical description of the dynamic motion of a 
ship requires six degrees freedom, three equations of force 
and three moment equations, these six equations are too 
complex for analytical solution. If only the rolling of 
the ship is considered, it may be modeled as a pure rota- 
tional motion about the longitudinal axis of ship, with no 
moments about the Y and Z axes. 
The ship in rolling, is analogous to a pendulum. The 
dynamic equation of a pendulum is: 
по+воч+ тын ф ЕЕ 
where m - mass of the pendulum. 
B - damping coefficient. 
g = gravity of the earth. 
1 = the length of the pendulum. 
$ = the instantaneous angle of pendulum. 
F = the excite force. 
Including hydrodynamic damping, the dynamic rolling 
equation of a ship is the same as the pendulum. 


аа т 


where J. = total moment of inertia around axis of roll. 
В. - hydrodynamic damping coefficient. 
K - static righting moment coefficient. 


ү = angle of drift, radians. 
KY = drift-angle coupling coefficient. 


ф = instantaneous angle of rolling. 
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The equation can be written as: 


T Ze 
ф % 262 9: % % шеф - Ш/ү 








| В 
where Es, 7 damping factor of the ship. -h 2 


E 
и 


$ `5 
natural frequency of the ship roll. x 
S 


ON 
О, = JT. 

When the ship is operating in a sea state, there is a 
multitude of forcing frequencies present and the integrated 
response to all these frequencies must be considered. Further 
experiments are therefore carried out with the model at zero 
speed in irregular beam waves. The rolling of the ship is 
caused by the force due to the wave slope. Then the equation 
becomes: 

P+ 2c. wo + флогф 


where w= the magnitude of wave slope. 


C. MODEL OF A SHIP 
The model of ship used for the example, to simulate and 
design the tank stabilizer is one given in Reference [6], the 


data and characteristics of the ship were the following: 


Displacement weight w = 936 tons. 
Metacentric height GM = 0.73 meter. 

Ship moment of inertia Je = 1290 tons-m-sec* 
Hydrodynamic damping coefficient В. - 60 tons-m-sec. 
Static righting moment coefficient Ko = 684 tons-m. 
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The mass moment of inertia of the tank stabilizer (J,) 
increased the value of 9 but it is very small and can be 
considered negligible when compared to Je > for convenience 
use the value of J. as above. 


The beam of the ship can be found by, 





2 2туСМ _ 2nVGM 
В = = EE 
та 5 
У. 


where C = a coefficient dependent upon the moment of roll 
and approximate is 0.797 — И 

Therefore the beam of the ship can be calcullated and the 

result is: 


B = 9.25 meters. 


D. PASSIVE TANK STABILIZATION 
ls Analysis 

A passive tank can only be designed for optimum 
performance for one ship condition of displacement and meta- 
centric height. When the ship condition varies, the free 
surface correction yu varies, together with the quantity of 
water required. It is therefore necessary to check the 
performance of the stabilizer under differemt ship loading 
conditions and it could be that a compromise design has to be 
adopted which is acceptable for all conditions. 

The free surface loss is generally expressed as: 


i 5 м В 


V 48 Ah 
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where i = the second moment of inertia of the free surface. 
У = the volume of displacement of the ship. 


w = the weight of active fluid. 


A = displacement of the ship. 
B = beam of the ship. | 
h = depth of fluid in tank. 

and Е - usual notation for free surface correction. 


The tank size parameter is limited primarily by 
static stability requirements. However, the tank should be 
designed with as large a free surface as possible consistent 
with the requirements of static stability. Choosing values 
of u is equivalent to varying the tank size. 

Since the pendulum мн is analogous to ship 
rolling, the effect on the ship of a hydrodynamic tank 
Stabilizer is analogous to the addition of a second pendulum. 
The combination of two pendulums in series, may be called a 


double-pendulum system as shown in Figure 5.1. 





Figure 5.1 Double-Pendulum System. 
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The first pendulum represents the ship rolling 
System, and the second pendulum represents the dynamic tank 
system. The length of the pendulums should be the same, so 
that the periods are nearly equal which means that both 
natural frequencies of the ship rolling and the tank are 
approximately the same. When the ship mass has been swung, 
the free oscillations are observed. When the tank mass has 
been attached, the masses of both ship and tank can be swung 
and the interplay and exchange of energy between the tank 
and the ship can be seen. 

» A The tank stabilization system uses water or some 
fluid mass device. There are many types of tanks. They can 
be designed as Doughnut tanks, completely filled sea ducted 
tanks, partially filled sea ducted tanks, U-tube tanks and 
free surface tanks. These tanks systems were suggested by 
Dr. J. H. Chadwick Jr. Typical tanks such as free surface 
and U-tube tanks have been suggested by G. J. Goodrich and 
Williams C. Webster respectively. Both types are described 
and designed in this chapter. 


A block diagram of ship and tank system is: 


ELS SHIP ? 


ТАМК CONTROL 





For a passive tank system the "Control block" 
does not exist, or may be considered a unit gain. Thus the 


block diagram for a passive system becomes: 


E SHIP : 
0 
| TANK | 


The relationship of the ship with a tank stabilizer 


can be represented as a double pendulum. The equations are: 


72% + B. + К. $ + 4.60 + Ке + F1 


Jono Ko + J,0 + B,O + K,o = 0 


where de = the moment of inertia of the ship rolling. 

В, - hydrodynamic damping coefficient of ship. 

К, - the static righting moment coefficient of 
ship. 

de = the moment of inertia of the tank. 

В, - the damping coefficient of the tank. 

К, - the static equivalent coefficient of the 
tank. 

dor - the coupling moment of inertia of ship/tank. 

Koy = the coupling static coefficient of ship/tank. 

F} = the force insert in ship as the wave slope. 
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2: Free Surface Tank Design 


The free surface tank as Figure 5.2 shown, is 
advantageous in decreasing the magnitude gain of rolling, 
and also decreases stability of the ship since the free 


Surface loss is high. 





Figure 5.2 Free Surface Tank. 


The equations of the ship and tank can be rewritten 


as: for the ship 
* : 2 4 2 20227772 


for the tank 


d 2 ec 2 T ° 2 А 
д 0.% - Ws > ин ox 25,0,9 + ш „0 = 0 


E 


where o, - natural frequency of the ship roll. 


E 
и 


t natural frequency of the tank. - 


aS 
ct jc u 


уч 
и 


damping factor of the ship roll. 
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- damping factor of the tank. 
g = gravity of the earth = 9.8 m/sec.* 
d = the distance between the water level in the 


tank and the roll center. 


V = wave slope function. 

u = the free surface correction factor. 
. 2 
1 _ 1 Ь 


V GM  2А6М. 

Both equations have the negative sign in the coupling 
terms since the restoring forces are in opposite directions. 
If the distance d is positive which means the level of water 
in the tank is higher than the rolling center, the magnitude 
gain is less as shown in Figure 5.3, so for convenience set 
d equal zero. 

Generally, the period of rolling of the tank (or 
natural frequency) must be designed nearly equal to the 


period of the rolling of the ship. The simplified equations 


become: 
$ + o. d ar od - 1 2,26 - wv 
6 + 0,0 + uo - 57 = 0 
where де = 250. 
a, = 25,0, 
ш = We = Wy 


By manipulation and transform, the system trans- 
mission function of the rolling angle due to the wave slope 


force is: 
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2 даа 2 
wp (S to. Stu, ) 


E. 3 2 2 
S %(а,4о,)5 + (дю аа, )5 


55) 2 274 
tu (a,ta,)Stw, (1-1) 

The solution of the equation can be expressed as 

an amplification factor A(w) and & the phase angle lag of 


the forcing function. 





42 
(ш) = lu 
At t y? 
А -1 Ү 
€ = tan X 
2 2 Е on (ur - ut] 
where X - (ші - ш - TONGE a Ж. 
n E 
ua, un 
Хо = ula, + ad 
(o^ -u*) tau 


These solutions can be used to investigate the 


effect of variation of the tank parameters a, and и. Calcu- 


t 
lations can be carried out using variation of each parameter 
in turn, keeping the other constant. Adjust both parameters 
to obtain the minimum magnitude gain. These values of para- 
meters are used to design the tank stabilizer. 
The computer simulation program used is due to 
Glavis [7]. The data of the ship in the previous section 
were substituted into the equation of ship and tank stabilizer, 


two unknown parameters inthe equation are o, and y as variables. 


tt 
Recall 
асы? 2 
$ : w (S ta Stw) 
еще NS a 22 2 
S *(o ta, )5 +(20,+0,0,.)5 +w, (a¿+a.)S+w, (1-1) 
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Numbering the data of the ship and parameters the 


equation becomes: 


de. 0.5302352+0.530230,5+0.28115 
54+(0.0465+а, )52+(1.06047+0.04650, )52+(0.02466+ 


0.53023, )5+0.28115(1-џ) 


For the sinusoidal signal the equation becomes: 


0.28115-0.530235^*j0.53023a,u 


T(jw) = ИН е ee ge 2 
ш -(1.06047*0.04651a,) 0 +0.28115(1-1)+j[(0.024662 


*0.53023a, )u-(0.0465*a, )u?] 


Then 


2 
M d k BM — E —cER- ME E 
[0^ - (1.06047*0.04651a,)u^«0.28115(1-u) ]^*[(0.024662 


(0.28115-0.53023w")“+(0.530230,1)* 


+0.530230,)w-(0.0465+a,)w"] 


Since the size of the tank depends on the free 
surface coefficient (u), if u is large, the tank will be big 
and the stability of the ship will be decreased because the 
metacentric height decreases. So the magnitude gain of 
rolling in the region of the natural frequency decreases but 
at low frequency the magnitude gain increases. 

Consider the computer results. Figure 5.4 shows 
curves of the tank damp ina coefficient versus frequency, with 
magnitude as curve index and the free surface coefficient is 
constant. On Figure 5.4 the free surface coefficient is 0.05. 
If a (tank damping coefficient factor) is approximately 0.25 


the height of the resonance peak is 13 db, which is the 
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smallest value obtainable for these conditions. This occurs 
at a frequency of about 0.72 rad/sec. At a frequency of 
about 0.18 rad/sec, the magnitude is only 1 db, and is even 
smaller at low frequency. 

On Figure 5.5, the free surface coefficient is 0.1, 
the line of constant tank damping coefficient factor is 0.35 
to give a resonance peak of about 11 db.which is the minimum 
obtainable. Note, however, that at lower frequencies the 
magnitude is higher than for the condition of Figure 5.4. 

If the free surface is set equal to 0.2 as on 
Figure 5.6, the tank damping coefficient should be 0.45, 
which gives 9 db for the resonance peak and 2 db at low 
frequency (0.5 rad/sec). Thus, if the free surface coeffi- 
cient is increased, the magnitude in the region of low 
frequency also is increased, though the height of the reso- 
nance peak is decreased. 

[he curves of Figure 5.7, ЭС 5.8 and Figure 5.9 
give tank damping coefficient factor versus magnitude. 
Frequency is the curve index, and the free surface coefficients 
are 0.05, 0.1 and 0.2 respectively. These Figures correspond 
to Figure 5.4, Figure 5.5 and Figure 5.6. By inspection they 
predict the same optimum values for the damping coefficient 
factor. 

Figure 5.10, Figure 5.11 and Figure 5.12 give the 
Bode diagrams (plot of the magnitude versus frequency) with 


the tank damping coefficient factor as curve index and the 
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free surface coefficient factor is constant. In this case 
the Bode diagram cannot show the best value of the tank 
damping coefficient factor unless that value has been plotted. 
However the Bode diagram is good for consideration of the 
values of the free surface coefficient for which the curves 
show different magnitudes at all frequencies. 

By inspection of the curve families the best values 
of the tank damping coefficient factor and the free surface 


coefficient factor are: 


the tank damping coefficient factor - 0.25 

the free surface coefficient factor - 0.05 
and/or 

the tank damping coefficient factor - 0.35 

the free surface coefficient factor = 0.1 


Note that the magnitude gain value for low fre- 
quencies increases with increases of the free surface (y), 
which means increasing tank size. This results from the 
losses in roll stiffness due to the tank free surface effect. 
It is also obvious that the roll response at low frequencies 
increases as tank size increases. 

For the free surface tank (as reference [8]), the 


relation between the parameters and the tank size is: 


в на 
Но ZAGM _ 
where ] = the length of the tank. 
b = the width of the tank can be equal or less than 
the beam of the ship. 
A = the displacement of the ship in tons. 
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GM = metacentric height 


The other suggestion equation is: 





42 2 
в = 2 
п 9 
where h = the depth of the water in tank. 


If the width of the tank is assumed to be equal to 
the beam of ship, the depth of the water in the tank is equal 
to 0.5 meter. If the selected free surface coefficient factor 
ms 0.1. 

1 = 0.82 meter. 

If the free surface coefficient factor 1s equal to 
0.05, the length of the tank is reduced to: | 

1 = 0.41 meter 
which can save space in the ship. And the other parameters 
such as Јр» В.» апа К, can be calculated. 

The time response to an initial condition is shown 
in Figure 5.13, for u » 0.0, 0.05, 0.1, and 0.2, and a, = 01205 
0.25, 0.35 and 0.45 respectively. The curve for u = 0.0, 
without antiroll tank stabilizer is very oscillatory, and is 
damped only by the damping of the ship itself. The other 
curves include the tank stabilizer, so the damping is higher, 
sand the curves reach steady state sooner. The time constant 
of the undamped system is about 42.5 seconds. With stabilizer 


added the setting time can be halved. 
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Figure 5.14 shows the time response to a sinusoidal 
input. The parameters have been varied as for Figure 5.13. 
The frequency of the input is set at about the natural fre- 
quency (0.728 rad/sec). For the ship without tank stabilizer 
the magnification is seen to be excessively large. When the 
tanks are added the system response is seen to be greatly 
reduced in magnitude. 

з. U-Tube Tank Design 

The U-Tude tank as designed in this chapter is 
assumed to be a rectangular tank. The mutual inertia term, 
Jet» which has been set to zero in preceding sections, is not 
necessarily zero. The value of J,, depends on the accelera- 
tion of the water relative to the tank, and this depends on 
the location of the tank. If the cross-duct passes below 


the center of rotation, y is a positive number which becomes 


st 
smaller as the tank system is placed higher, as shown in 
Figure 5.15a and b. 


Consider that the two parameters J, and J,, are 


given by: 
2 5 
J, 7 pAg R $ Ap ds 
0 А. 
2 5 
Jog = PAg R ) 0, 45 
R 
where p = mass density of the fluid in tank. 


- 


Ag - cross-sectional area of vertical leg of the tank. 


КЁ = lever arm of the tank. 
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р. = perpendicular distance from "virtual" center of 
rotation to the tangent at any point on the 
trajectary. 

A_ = the cross-sectional area at any point 5. 


S 
These two equations can be related as: 





dot _ 9 р, А, 
йы Б pz i ads 
t 0 0 
2D A 
- 2Һ 4 —[1—— 
0 


and h = the height of the water in the tank. 
Note that р, can be negative if the tank is placed 


above the center of rotation. So Jet can be zero if: 





h A 
a 0 
-D, — A 
S 
Recall | 
9 9 3 В; $ = К. $ P 9569 P Kot? 4 Ko sy 
J;u6 * K ¿0 + J¿0 + B,0 + ке - 0 


By using the definition as before the transformed 
equations become. 


J 


2 2 51229722 t _ 2 
(S ta, Sto.) + (- тас —)0 = шеф 
5 S 
( - “st 52-09) + (544 S+ ar: = 0 
J, n aj An. 
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К 
Define үес В (зи as in the free surface tank) 
S 
J 2 пад 
75" - Үй- 2% + ~ 5 
t 0 
a, во 
then "st = BX 
J 2 
5 ши 


Using Cramer's rule, the ship roll transfer function 


due to the wave slope force is: 


we (S240 5409) 

£(s) = ; n man 
Х В, „4 2 2 2 2 4 2 
Due +(a,ta,)S +(20,+0,0.-2X8)S tw Carta) Ston Bu, 


Inserting the values of the parameters of the ship, 


the equation becomes: 


0.5302352+0.5302305+0.28115 


als) J ПЕ 1 3 2 
[(1-1.886X°B)S '+(0.0465+0)S~+(1.060474+0.0465-2X8)S“+ 
. (0.02466-0.530230)540.28115-0.530238 


then 


ES (0.28115-0.53023w%)%+(0.530230w)* 


[(1-1.886X^8)u* -(1.0604740.0465-2X8)0^*0.28115-0.530238]^ 
+Г(0.024662+0.53023а)0-(0.0465+а) a7 ]^ 


Figure 5.16a, b and c are plots for X-3 and 820.05, 
there are two resonance peaks in the magnitude curve, one 
for the ship itself and the other for the tank. The value of 
a at the minimum peak magnitude is about 0.45 and only one 
peak magnitude is obtained. The point P in Figure 5.16c 
corresponds to the line of -2 db in Figure 5.16a. 
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If the tank is placed as high as possible, the value 
of X may be zero or negative. If X=0, the result is the same 
as for the free surface tank with es. The plot of X=-1.0 
and 820.05 is shown as Figure 5.17a, b, and c. The points P 
and Q on Figure 5.17c are the intersection of the a lines, 
these points correspond to the lines of the frequencies 0.6 
and 0.9 rad/sec on Figure 5.17a. In this frequency region, 
if the value of a is decreased, the magnitude also decreases. 
Figure 5.17b shows that minimum peak value of magnitude is 
obtained with a=0.7. 

The U-Tube tank design equation is: 

К. = 2рдА К 
where P = mass density of the water. 

g = gravity of the earth. 
А, - cross-sectional area of each vertical leg of the 
tank. 
R = distance from ship center line to center of a tank 
Vertical leg. 

The product of p and g is unity for the water, so both 

can be eliminated from the equation. 


And the cross-sectional area of the channel between 


tanks is given by: 


2 

_ Bo Ко С 
АР E 2 
д - hw. 
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Frequency Parameter Plane (Bode Diagram). 
free-surface tank. 


Figure 5.3 
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Figure 5.4 Frequency Parameter Plane for Free-surface Tank. 
Plot of Tank Damping Factor vs Frequency with 
Free-surface, u-0.05 and M as curve index. 
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Figure 5.5 Frequency Parameter Plane for Free-surface Tank. 
Plot of Tank Damping Factor vs Frequency with 
Free-surface, u-0.1 and M as curve index. 
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Figure 5.6 Frequency Parameter Plane for Free-surface Tank. 
Plot of Tank Damping Factor vs Frequency with 
Free-surface, u-0.2 and M as curve index. 
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Figure 5.7 Frequency Parameter Plane for Free-surface Tank. 
Plot of Tank Damping Factor vs Magnitude with 
Free-surface, u-0.05 and o as curve index. 
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Figure 5.8 Frequency Parameter Plane for Free-surface Tank. 
Plot of Tank Damping Factor vs Magnitude with 
Free-surface, џ=0.1 and o as curve index. 
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Figure 5.9 Frequency Parameter Plane for Free<surface Tank. 
Plot of Tank Damping Factor vs Magnitude with 
Free Surface, y = 0.2 and w as curve index. 
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Frequency Parameter Plane for Free-surface Tank. 
Bode Diagram, Plot of Tank Damping Factor as 


curve Index with Free-Surface, u 


Figure 5.10 
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Frequency Parameter Plane for Free-surface Tank. 
Bode Diagram, Plot of Tank Damping Factor as 
Curve Index and Free- СВ, u=0.1, 
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Figure 5.12 Frequency Parameter Plane for Free-surface Tank. 


Bode Diagram, Plot of Tank Damping Factor as 
Curve Index and Free-surface, и-0.2. 








Figure 5.13 Initial Condition Response with Different Value 
of Tank Damping Factor and Free-surface. 
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Figure 5.15 U-Tube Tank 
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Figure 5.16a Frequency Parameter Plane for U-Tube Tank. 
Plot of Tank Damping Factor vs Frequency 
with X = 3, 8=0.05 and M as Curve Index. 
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Figure 5.16c Frequency Parameter Plane for U-Tube Tank. 
Bode Diagram, Plot of Tank Damping Factor 
as Curve Index with X = 3 and 8720.05. 
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Figure 5.17a Frequency Parameter Plane for U-Tube Tank. 
Plot of Tank Damping Factor vs Frequency 
with X = -1, ß=0.05 and M as curve index. 
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Figure 5.17c Frequency Parameter Plane for U-Tube Tank. 
Bode Diagram, Plot of Tank Damping Factor 
as Curve Index with X » -1 and 870.05. 
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This is for the rectangular U-Tube tanks. With the 
equations as above, the values of cross-sectional area of 
the vertical leg of the tank, АЛТ the cross-sectional area 
of the channel, A.» can be calculated and tank designed. 

The depth of the water in tank, h, and the distance 
between the rolling center and the channel of the tank, Do > 
can be calculated by: 

2 р, A. 


X=2h + Д 


0 

In the case of X = 3, if the tank is located below 
rolling center, h should be less than 1.5. 

E.g., if X = 3, B=0.05, the tank damping factor is 


0.45. Then the tank size is: 


if R = 4, the width of the tank is about 1 meter. Therefore 


128222 013 15 


and if h = 1.3 meters, the result of pe 1s 0.4686 meter“ 


and р, - 0.859 meter. 


БА ACTIVE TANK STABILIZATION 
17 Analysis 
Activated-antiroll tanks are considered in detail. 
Particular attention is paid to the design of the tank control 
system. The equations of motion are the same as those of the 
ship with passive tank stabilizer, but must be modified to 
include the pump used in the activated tank. Active tank 


control can best be selected on the basis of minimizing the 
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рт 


response of the ship to an impulsive roll velocity, and it 
will be shown that the active tank is significantly better 
in reducing the roll of ship than the ordinary passive antiroll 
tank. The active tank stabilizer studied uses the U-Tube tank 
design. | 

Active tank analysis parallels that of passive tanks 
and/or the double pendulum. The tank system is activated by 
inserting a force (such as a pump system) into the second 
equation of passive tank system. The equation can then be 


written as: 


ship: 72% + В. $ + К.ф t J t8 4 К. „0 = к vy 


! 
“TI 


tank: 9560 i Kee? + J,6 + Ч, t K, 6 р 


where Fp = the force due to pump system. 
The effect of the roll sensing control on the force 
produced by the pump system can be expressed as: 


$+6.5°)% 


К (0,540 
_ |2222 3 
F5 s) Е Ши, 


where к ET gain of the control pump. 
б - roll angle sensitivity. 
б. - roll rate sensitivity. 
G3 - roll acceleration sensitivity. 
t = time delay constant. 
For convenience, we set т=0 in the computation, since 


the time delay is small. 
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The active tank system can be shown in block diagram 


form as: 


PASSIVE 


20% ф 


РЫМР 


ЯВ Po | | Кр [CONTROL 


where Е. (5) - transfer function of the ship rolling due to 







wave slope alone which is the same as for the 
passive tank system. 
(5) - transfer function of the ship rolling due to 


pump alone, which can be expressed as: 


А В A 
B. MENO 3.,,2 2 2 4 2 
D +(а +9, )5 *(2u, ta, a, -2X8)S tw (а, +9. ) 5+0 „Вр 
5 
292 
= z5 


2. System Design 
From the analysis in this section, the system per- 
formance can be studied using computer simulation. The pump 
coefficient factor к, (к, = Ту) is one variable parameter 
which may be called а. It is varied as the pitch ratio of 
the propeller of the pump, and can be assumed to be varied 


from zero to one and zero to minus one. For convenience the 


actual system may be designed with only positive variation and 
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the direction of the water motion between tanks is controlled 


by valves. 


From the block diagram the system transmission 


function is: 





ш2[52+а $4ш2 1 27 y 
S(s) mos 5 n t n цэн ЕЕЗ: | 
XT о MM Me 2 2 > 
> +(a,ta,)S +(Zw ta,a, 2XB K,638)5 + (ша, ша, 
n 
1 4 2 ! 
-К,6)5%ш,-Вш,-К,618 
" 2 
where F (5) = K y (67 +6,5+635 )ф. 
Define X - 3 
B - 0.05 
ын = 0.4 
1 
and Ko = а 


These values for a passive tank gave a peak resonance 
of minimum magnitude at about 10 db at 0.6 rad/sec and is 1 db 
at the lower frequencies. Using these values the equation 
becomes: 

(s) - с 0.5302357+0.2120925+0.28115 
0.15135 +0. 44655 +(0.7791-0.0506,)S %(0.237-0.0506.)5 


+0.25464-0.0506, 


Figure 5.18 shows the Bode diagram for 6721, 
G,=6,=0. Ін а 15 more positive, the frequency and the magni- 
tude of the resonance peak are slightly decreased. At the 
low frequencies, the magnitude is increased. If a is more 


negative, the result is reversed. To obtain the minimum 
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Figure 5.18 Frequency Parameter Plane for Activated Tank. 
Bode Diagram, Plot of Pump Coefficient Factor 


as Curve Index (6,71, 6576.70) and X23. 
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magnitude gain, at the lower frequencies from 0 to 0.6 rad/sec 
а 15 -1, for the higher frequencies a must be +1. So, the 
pitch ratio of the propeller of the pump may be kep constant 
and the direction of the water motion reversed by switching 
the valves. Then the maximum magnitude is 7 db at 0.6 rad/sec 
and -1 db at the lower frequencies. 

The characteristics of the active device (pump, 
valves, sensor) can be treated as the variation of the value 
of the parameter a versus frequency. This value is required 


to switch at the resonant frequency as shown in Figure 5.19. 






w rad/sec 


-] 


Figure 5.19 The variation of control device. 


(X23, @1=1, 6576470) 


For G,=0.5, 6,=1.25 and 6,=2.4 the results are shown 


1 2 3 

as Figure 5.20a, b and c. Figure 5.20a is a plot of a versus 
frequency, and magnitude is the curve index. When a is more 
negative the magnitude is seen to be smaller. If a=] the 


maximum peak magnitude is obtained at the frequency of about 
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Figure 5.20a Frequency Parameter Plane for Activated Tank. 
Plot of Pump Coefficient Factor vs Frequency 
with M as Curve Index (6170.5, G571.25, 
6472.4) and X=3. 
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Figure 5.20b Frequency Parameter Plane for Activated Tank. 
Plot of Pump Coefficient Factor vs Magnitude 
with w as Curve Index (6120. 5, Go=1.25, 


6.=2.4) and X=3. 2 
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Figure 5.20c Frequency Parameter Plane for Activated Tank. 
Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (6170.5, 651.25, G3=2.4) and 
X=3. 
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0.6 rad/sec. These results correspond to Figure 5.20b which 
shows the frequency curve. At the same frequency of 0.6 
rad/sec, both figures show maximum magnitude for all values 
of a from -1 to +1. Figure 5.20c is the Bode diagram and 
shows that there is no phase lag between ship and active tank, 
and if a is more negative the magnitude is smaller. So, at 
а=-1, it produces the smallest magnitude at all frequencies, 
with the peak equal 5 db at 0.6 rad/sec. 

This result seems to be the best since all the 
parameters of the pump system (pitch, sensing gain control, 
and valves) are fixed. 

From both results as mentioned above, the system 
can be designed as a linear device. 

If there is no limit on the height at which the 
tank can be located, it is possible to choose the height, so 
that Х-0. If this is done, by inspection the value of 8720.05 
and a,=0.3 are obtained. Substitute these values into the 
equation of the system which becomes: 

2 


(5) O20 30235020 259069550425 ho 


$1+0.34655°+(1.07436-0.05К_С (5°+(0.18373-0.05K ,6,)S 


p 3 


*0.25465-0.05& G4 


Figure 5.21a and b show the results with 6771 and 
G 


=6.=0. Figure 5.21a is a plot with magnitude as curve 


e 3 
index; point Q is the maximum magnitude with a=-0.65 at a 
frequency of 0.728 rad/sec. If a=0.15, the maximum magnitude 


of 9 db is obtained which is the minimum peak. Figure 5.21b 
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Figure 5.21a Frequency Parameter Plane for Activated Tank. 
Plot of Pump Coefficient Factor vs Frequency 


with M as Curve Index (6471, 6576470) and X-0. 
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Figure 5.21b Frequency Parameter Plane for Activated Tank. 
Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (6.71, 6:6: =0) and X=0. 
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shows the Bode diagram. To obtain the minimum magnitude, in 
the frequency region between 0.6 and 0.85 rad/sec the value 
of a must be maintained at *1 and for the other frequencies, 
at -1. There are two magnitude peaks of approximately 7 db 
at the intersection of а=+1 апа а=-1. Рог the lower frequen- 
cies the magnitude is -1 db. 

From the above result the active device parameter 
value, a, can be plotted versus frequency as in Figure 5.22. 
This figure shows two values of a at the frequencies between 
0.6 and 0.85, the straight line corresponds to the minimum 
magnitude, and the curved line corresponds to a constant 


magnitude in that region. 
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Figure 5.22 The Variation of Active Device 


(X=0, G,=1, G,=6,=0) 


Figure 5.23 shows the Bode diagram of the system 
for G,=0.5, Go=1.25, and G,=2.4. For minimum magnitude with 
the control device, the required variation in a is shown in 


Figure 5.24. 
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Figure 5.23 Frequency Parameter Plane for Activated Tank. 
Bode Diagram, Plot of Pump Coefficient Factor 
as Curve Index (6170.5, 6571.25, 6472.4) and 
Х-0. 
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Figure 5.24 The Variation of Control Device 
(X=0, G,-0.5, G,-1.25 6.52.4) 
Figure 5.25, which is the Bode diagram, compares 

the results for an undamped ship and for a ship damped by a 
tank. All have only one resonance peak. The magnitude of 
roll in the undamped ship at resonance is the highest, damping 
clearly reduces the amplitude. However, in the passive tank 
system the magnitude at lower frequencies is slightly higher, 
and the curves show that the U-Tube tank is more advantageous 
than the а е The active tank system gave the 


most advantageous behavior as shown. 
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Figure 5.25 Compare the Results in Magnitude Gain. 
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VI. CONCLUSION 


Analysis techniques can be employed effectively in system 
design particularly when an exact analytical solution is not 
available and where a trial and error method is the only 
solution possible. 

The frequency parameter plane techniques have been used 
in this thesis to analyze and design a system. The system 
elements which may be linear or nonlinear were represented 
by adjustable parameters. То select a stable system, the root 
parameter plane was used and further mapped onto the frequency 
parameter plane. The goal was to obtain the best set of para- 
meters for a stable system. 

Rather than adjust parameters for every trial, one graph 
was made and the parameter value corresponding to the best 
curve was chosen. This appreciably reduced the number of 
trials and the computation time in the overall design. 

In the system which is presented the tank stabilizer 
design might not be optimal, however further effort on para- 
meter adjustments by iteration of the described technique 
would improve the design. | 

The tank stabilizers of a ship were designed in this 
thesis by Tos frequency parameter plane technique, they 
were very successful and effective. If the relationship 
between the designed parameters and tank size are true, the 


result of designing by using the method would be feasible and 
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believable. In passive tank design, the tank size depends 
on the free surface factor. The tank damping factor depends 
on the skin of the material of the tank, which would be tested 
to obtain its value from the analysis. The value of X in 
U-Tube tank would have to be chosen within the range which 
maintains ship stability. 

In the activated tank design, the gain control constants 
were varied by trial and error to adjust the damping.  Phys- 
ically only one parameter can be varied, the pitch radio of 


the propeller of the pump. 
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